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We propose an experiment for entangling two spatially sep- 
arated Bose-Einstein condensates by Bragg scattering of light. 
When Bragg scattering in two condensates is stimulated by a 
common probe, the resulting quasiparticles in the two conden- 
sates get entangled due to quantum communication between 
the condensates via probe beam. The entanglement is shown 
to be significant and occurs in both number and quadrature 
phase variables. We present two methods of detecting the 
generated entanglement. 

PACS numbers: 03.75.Fi,03.65.Ud,42.50.Dv 



Inseparability of quantum states of two or more sub- 
systems is the most significant feature of quantum me- 
chanics. Apparently puzzling, yet most profound, first 
formulated as a paradox Qj, this inseparability known 
as quantum entanglement lies at the very heart of non- 
classical physics. Further, as a basic resource for quan- 
tum information processing, it has become a focal theme 
of research in modern physics and many issues in the 
foundations of quantum mechanics. Generation and ma- 
nipulation of entanglement is, therefore, of prime inter- 
est. Bose-Einstein condensates (BEC) Q of weakly in- 
teracting atomic gases seem to be suitable macroscopic 
objects for producing many-particle entanglement A 
BEC has intrinsic entanglement character due to reduced 
quantum fluctuations in momentum space. For instance, 
in the condensate ground state, a pair of mutually oppo- 
site momentum modes is maximally entangled in atomic 
number variables 

Stimuated resonant Bragg scattering of light by a con- 
densate generates quasiparticles ||, predominantly in 
two momentum side-modes q and — q, where q is the 
momentum transfered from light fields to the atoms. Mo- 
mentum side- modes are the excited states of a BEC, 
atoms in such a state collectively behave as quasipar- 
ticles. Bragg spectroscopy || with coherent or classical 
light produces coherent states of the quasiparticles in a 
BEC. When these quasiparticles are projected into par- 
ticle domain, that is, into the Bogoliubov-transformcd 
momentum modes 0, they form two-mode squeezed as 
well as entangled state Q] . Bragg spectroscopy with non- 
classical light can generate tripartite entanlement Q in 
a condensate. In addition to atom number and phase 
variables, spin degree-of-freedom of a spinor BEC § can 
be useful in describing entanglement in spin variables. 
Thus, BECs offer a fertile ground for studying different 
aspects of entanglement. Apart from BECs, multi-atom 



entanglement in other macroscopic systems has been re- 
alized H on the basis of collective spin squeezing jl0|,[ll| ■ 
Further the entanglement in collective spin variables of 
two ensembles of gaseous Cs atoms has been experimen- 
tally demonstrated ]l2] . Continuous variables like the 
quadratures of a field mode (which are analogous to po- 
sition and momentum) have also been employed ]l3[ ] in 
entanglement studies. 

We here propose a scheme for producing quantum en- 
tanglement between two spatially separated BECs of a 
weakly interacting atomic gas. The entanglement we con- 
sider is in quasiparticles of BECs. The proposed exper- 
iment is schematically shown in Fig.l. The condensates 
A and B are illuminated by pump lasers LI and L2, re- 
spectively. A single stimulating probe laser L3 passes 
through both the condensates. All these three lasers are 
detuned far off the resonance of an electronic excited 
state of the atoms. The frequencies and the directions 
of propagation of these lasers are so chosen such that 
Bragg resonance (phase matching) conditions of scatter- 
ing in both the condensates are fulfilled. The Hamilto- 
nian of the system is H — Ha + Hb + Hp + Haf + He f ■ 
Retaining the dominant momentum side-modes q and 
-q only under Bragg resonance condition, in the Bogoli- 

ubov approximation jjj, Ha — tiuj^ ^a q <S q + aL q a- q ^ , 
where d q represents quasi-particle with mometum q, and 
ujB = [(uj q + g ) 2 - (g ) 2 ] 1/2 is the frequency of Bogoli- 



ubov's quasi-particle 0. Here 



h A q l 
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J 1 2m ' f Ira 

is the chemical potential with £ = (87rnoa s ) -1 / 2 being the 
healing length. Similarly, H B = hto^ (/3 q /3 q + /3L q /3_ 




FIG. 1. The scheme for creation of entanglement. A and B 
are two condensates, LI and L2 are pump lasers, L3 is a com- 
mon entangling probe laser. Both the pumps have same wave 
vector ki, probe's wave vector is It2. The probe is red-detuned 
from the pumps. The lasers are in Bragg resonance with a 
particular momentum mode q of both the condensates. 

wih /3„ being the quasi-particle operator of the con- 
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densate B. The pumps are treated classically. Let c 
represents the common probe field mode, then the field 
Hamiltonian Hp — —TiS&c, where S = u>\ — u>2 is 
the pump-probe detuning. The quasiparticle operators 
&($) are related to the particle operators a(b) by Bo- 
goliubov's transformation : a q = u q ctq — v q 6^_ v where 
v q = {u 2 q - 1)V2 = - 1)] 1/2 . The atom-field 

ineteraction Hamiltonian for condensate A is 

Haf = hr]AC ! (al l + a-q) + H.c. (1) 

where t\a = V^A^Afq is the effecive atom-field coupling 
constant. Here Na is the number of atoms in condensate 
A , Qa is the two-photon Rabi frequency of an atom in A 
and f q = u q — v q . Hbf is given by the similar expression 
as Haf with subscript A replaced by B and a replaced 
by (3. The Heisenberg equations of motion are 

d q = -iujqCtq - ir)AC f (2) 
a_ q = iWgdLq + in A c ] (3) 

c = -iStf + iivAi&q + al q ) + VB((3q + #L q )] (4) 

The Heisenberg equations of /3 q and /31 q are similar to 
those of d, but d and 77,4 should be replaced by and 
t)b, respectively. 

We next discuss how to quantify entanglement between 
two BECs. If the entanglement occurs in number oper- 
ators of the quasiparticle modes 1 and 2, then it can be 
quantified by the parameter ]lO|,pj 

£„(1, 2) = ([A(m - n 2 )] 2 )/(( ni ) + (n 2 ». (5) 

If < 1) then the two modes are entangled. If the en- 
tanglement is described by two noncommuting Gaussian 
operators X and P which are analogous to position and 
momentum variables, then the entanglement parameter 
is defined by Q 

e P a,2) = l -\mx 1+ x 2 )f) + ([a(p x - p 2 )] 2 )] (6) 

The two modes are entangled in quadrature phase, when 

For numerical illustration, we consider two homoge- 
neous identical Na condensates. We here enlist the im- 
portant results: 1) If the modes qi of A and q 2 of B are 
in Bragg-resonance with the the respective Bragg pulses, 
and if the effective coupling of B (773) is stronger than 
that of A, then entanglement arises between q x of A and 
— q 2 of B only, other pairs of modes are immune to any 
entanglement. In fig. 2, we display entanglement param- 
eters between these two chosen modes as a function of 
time. We set qi = q 2 = q. The effective coupling can be 




t(ns) 

FIG. 2. Entanglement parameters £ n and £ p between 
q-mode of A and -q-mode of B as a function of time in fis. For 
both the condensates, uif = 0.21 MHz, q = 2£ -1 and S = 0.17 
MHz. The coupling constants r\A = 1.62 MHz, r\B = 1.25r/i. 
Both the condensates are initially in the ground states (vac- 
uum of quasi-particle operators) , and the common probe field 
is in coherent state with average number of photons equal to 
10. 

made different either by using pump lights of different 
intensities or taking different atom numbers for the two 
otherwise identical condensates. 2) For equal couplings, 
there is no entanglement between any pair of modes. 
Fig. 3 shows the variation of entanglement parameters as 
a function of the ratio of the two coupling constants at 
a fixed time. 3) We find entanglement both in quasi- 
particle (phonon) modes (d,/3), and in the Bogoliubov 
transformed modes of quasiparticles which we call parti- 
cle or atomic modes (a, b). However, in atomic modes, 
entanglement is weaker than that in quasiparticle modes. 
It is worth mentioning that in a single condensate, as 
shown in Ref. coherent light scattering can generate 
entanglement only in atomic modes, and not in phonon 
modes. In contrast, one can generate entanglement in 
phonon modes in two condensates by coherent light scat- 
tering. The light scattering events occuring at A and B 
are not independent, since a quantum communication has 
been set between the generated quasiparticles in A and 
B via the common probe. Had we treated the common 
probe classically, then the Hamiltonian (Eq.(l)) would 
have been linear in the atomic operators. A Hamiltonian 
linear in Bosonic opeartors can not generate nonclassical 
correlation. Therefore, the probe must be treated quan- 
tum mechanically. The probe carries with it quantum 
fluctuations of one condensate and transfers a part of it 
to the other leading to the entanglement between the two 
condensates. 

To explain the results further, we here resort to an ap- 
proximate analysis. Let us suppose, cu q « t]a(b) an d 
5 « r]A(B)i then we can neglect the diagonal terms pro- 
portional to uj q and 5 in the Hamiltonian. For equal 
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FIG. 3. The entanglement parameters £ n (solid) and £ p 
(dashed) as a function of tib/va for a fixed time t = 0.75 /xs. 
The other parameters are the same as in Fig. 2. 

coupling (r/A = t/b), from Hcisenberg equation of mo- 
tions, it then follows that a q (t) + /3_ q (i) = a q (0) + 
/?l q (0), that is, the superposition operator S = a q + /3^ q 
becomes a constant of motion. Let us write the quadra- 
tures X A = -i=(d q + d q ), P A = -i^i&q - a q ), and 
similarly for Xb and Pb- Then it can be shown that 
£ P = i[([A(S + St)] 2 ) - ([A(E - St)] 2 )]. For equal cou- 
pling and the initial states being in vacuum or in coher- 
ent states, one obtains £ p = 1, that is, the two modes are 
unentangled. Let us then consider the case of different 
couplings, for short times characterized by rjAt « 1, and 
i]st « 1, we obtain perturbative solutions of a q (i) and 
(3-q(t) upto the second order in time. Using these solu- 
tions, we calculate £ p = 1 — T] A riBt 2 (i — which is less 
than unity (the two modes are entangled in quadrature 
variables) if ^a^b^I — ^) > which is only possi- 
ble if rjA 7^ rjB and t/a < t]b- Similarly, we can prove 
that for -q (off-resonant) of A and q (resonant) of B, 
£ p = 1 + (rjBt) 2 (1 — —) which is always greter than unity 
for t]a < r/B- For the same resonant q-mode of A and 
B, £ p = 1 + T] B t 2 /2 + {tf A + ry|) 2 i 4 /4, which is always 
greater than unity. In the same way, we can show that, 
for the remaining mode-pair (— q, — q), £ p is also greater 
than unity. 

Next, we prove that, to generate entanglement in num- 
ber variables (£ n ), the two coupling parameters should 
also be different. Substituting hi — a q a q and h 2 = 

/3l q /3_ q in Eq.(5) and using the pertubativc solutions we 
can express = 1 — R/{{h\) + (n 2 )) where 

A = 8»&t 4 fo| - 2»& + 4n p ( V 2 B - Va )} (7) 

where n p is the initial number of photons in the co- 
herent probe beam. Now, < 1 implies that R > 
which amounts to (?7b/?7a) 2 > 1 + 1/(1 + 4n p ), that is, 
Vb > Va- On the other hand, if rj B < f]A, £n > 1- We 



also carry out an aletrnative analysis to check whether 
the two resonant modes q of A and B exhibit any en- 
tanglement in other parameter regimes. By neglecting 
the off-resonant mode — q in both the condensates and 
keeping only the resonant mode, it can be analytically 
proved that £ p (q, q) = 1 + sinh 2 (r/t) and £ n (q, q) = 

1 + C 1 + ^) §j^t? sinh2 (^)' that is > both the Pa- 
rameters £ p and are always greater than unity. Here 

V = VVa + Vb- 

We next show how a set up as shown in Fig. 4 can be 
utilized to verify the generated entanglement. After the 
process of generation of entanglement, the duration of 
which can be typically on the order of 1 to 100 fis, is 
over, the lasers LI, L2 and L3 are switched off. Two 
different pairs of verifying pump-probe Bragg pulses are 
applied to the condensates, as described in the caption of 
Fig. 4. The two probes should be derived from a common 
source. The modes q of A and — q of B arc in Bragg 
resonance with the respective Bragg pulses. The effective 
field-condensate couplings for both the condensates are 
very small compared to the Bogoliubov frequency uj^ . 
Let c pro be,A and c pro b e ,B denote the verifying probe field 
modes for the condensates A and B, respectively. By 
neglecting the off-resonant terms ajq^q and a! 2q ci:_ q in 
the Hamiltonian, the time evolution of the output probe 
modes, in a frame rotating with pump-probe detuning <5, 
can be written as 




FIG. 4. The scheme for verification of the entanglement. 
Apart from pump lasers, two extremely weak verifying probes 
of the same frequency as that of entangling probe - one each 
for each condensate - are switched on . The probes acting 
on A and B have momentum k2 and — k2, respectively. The 
momentum of pump laser for A is ki, while that for B is — ki. 
Thus, the mode q and — q are in Bragg resonance with A 
and B, respectively. The output from the two phase sensitive 
detectors (PSD) can be integrated by an integrator followed 
by another phase-sensitive detection of the integrated signal. 

cp°robe b ^ s gi ven by similar expression as above with a 
replaced by (3; and the subscripts A and q replaced B and 
— q, respectively. Thus the output probes have oscillating 
parts at frequency 5 — uo^ proportional to the quasiparti- 
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cle amplitudes a q and /3- q , and at frequency S + uj^ pro- 
portional to the amplitues a_ q and /3 q . Therefore, phase- 
senitive measurements of the spectral components of the 
output probe beams corresponding to these frequencies 
would provide measures of the quasiparticle operators. 
The output from both the PSDs can be integrated by an 
integrator and the integrated signal can also be measured 
by another PSD. By repeating the same measurements 
under identical conditions we could calculate the number 
variances or correlation functions of interest, which can 
be employed to calculate the entanglement parameter in 
number operators, i.e., £„. For calculating entanglement 
parameters in quadrature phase variables, both the out- 
put probe beams coming from A and B, can be mixed 
via a beam splitter to form the superposition operators 
£ which can be measured by a similar phase sensitive 
detection scheme (not shown in Fig. 4). 

Following the recent experiment of Ketterle's group 
|l5|| , we also suggest that the quasiparticles can be de- 
tected by imparting a large momentum to them with ad- 
ditional Bragg pulses. Alternatively, in the large mo- 
mentum regime (q » the Bragg-scattered atoms 
which essentially behave as free particles (ujg oc q 2 ), can 
be outcoupled by swithcing off the trap. Since entan- 
glement is between two opposite momentum states, by 
proper geometric arrangement, the two moving entangled 
atomic ensembles can be made to collide and interfere. 
From the interference pattern obtainable via absorption 
imaging, the atomic number fluctuations can be deduced 
using the theoretical model used in Ref. jlfj], and thus 
entanglement parameter in number variables can be cal- 
culated. 

In conclusion, we have theoretically demonstrated how 
light scattering leads to quantum entanglement between 
two Bose-Einstein condensates. We find that the quasi- 
particle or phonon as well as free-particle momentum 
modes of two condensates can be entangled. The quasi- 
particle state can be sufficiently long-lived due to weak 
nature of interatomic interaction and the constraints im- 
posed by momentum conservations. The generated en- 
tanglement may be useful in quantum communication us- 
ing coherent light p^| . We have particularly focused on 
the conditions under which the entanglement can be ob- 
tained. We have also suggested how quasiparticles could 
be studied by using Bragg scattering of far off resonant 
fields. 
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